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Abstract 

In this paper, the phase diagrams and the critical behavior of the spin-1/2 anisotropic XXZ 
ferromagnetic model (the anisotropic parameter A € (— oo, 1]) on two kinds of diamond-type 
hierarchical (DH) lattices with fractal dimensions df = 2.58 and 3, respectively, are studied via the 
real-space renormalization group method. It is found that in the isotropic Heisenberg limit (A = 0), 
there exist finite temperature phase transitions for the two kinds of DH lattices above. The systems 
are also investigated in the range of — oo < A < and it is found that they exhibit XY-like fixed 
points. Meanwhile, the critical exponents of the above two systems are also calculated. The results 
show that for the lattice with df = 2.58, the value of the Ising critical exponent vi is the same as 
that of classical Ising model and the isotropic Heisenberg critical exponent is a finite value, and 
for the lattice with df = 3, the values of vi and z^h agree well with those obtained on the simple 
cubic lattice. We also discuss the quantum fiuctuation at all temperatures and find the fiuctuation 
of XY-like model is stronger than the anistropic Heisenberg model at the low-temperature region. 
By analyzing the fluctuation, we conclude that there will be remarkable effect of neglecting terms 
on the final results of the XY-like model. However, we can obtain approximate result at bigger 
temperatures and give qualitatively correct picture of the phase diagram. 
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I. INTRODUCTION 



The spin-1/2 anisotropic Heisenberg model, or the XXZ model, is one of the most im- 
portant quantum spin models, which can be used to explain the quantum essence of many 



magnetic materials, such as La2Cu04 and thin films of "'Hejll. Recently, for this model, 
many works have been done on translational symmetric latticesj^S]. For example, F C 
Alcaraz and A L Malvezzi studied the XXZ chain in the presence of an external field and ob- 
tained the phase diagram of this system accurately [3|. Continentino and Sousa, respectively. 



have investigated the critical properties of the two-dimensional (2D) anisotropic Heisenberg 
modelj4-^, and found that, in the isotropic Heisenberg limit, this system does not exhibit 
finite temperature phase transition (i.e., the critical temperature = 0), which is in accor- 
dance with the theorem of Mermin and Wagner 91]. However, as to this model on the simple 



cubic lattice Tq no longer equals zero in the isotropic 



10- 



ieisenberg limit, which means the 



12|. 



existence of finite temperature phase transition [7|, |8, 

On the other hand, there has been some interest in the critical phenomena of XXZ model 
on fractal lattices, and particularly on the diamond-type hierarchical (DH) lattices. Because 
of the special geometrical and topological property, DH lattices are good candidates to inves- 
tigate the spin systems in non-integer dimensions, and since they have a much lower symme- 
try than other fractals, so they may provide insights into other low-symmetry problems such 
as random magnets, surfaces, and the like 13|, |2J]. In these aspects, some theoretical works 
have been done. In 1983, using the renormalization group (RG) method, Caride et al have 
investigated the critical behavior of the Heisenberg model on the Wheatstone-bridge-basis 
hierarchical lattice. Their results show that in the isotropic Heisenberg limit Tc approaches 
zero as a continuous function of the anisotropic parameter ij] . Latterly, the phase diagram 



of this model was obtained by Souza on another kind of DH lattice with fractal dimension 



d 



f 



15|, and the results agree well with those on the square lattice. Recently, using a 



real-space RG method, the anisotropic Heisenberg spin-glass model on a three-dimensional 



DH lattice has been studied 



16|- 



In the study of the anisotropic Heisenberg model, many effective approximate methods 



have been a 
method 



m 



ied. Such as the mean-field approximat ion 1 1 0| . series expansion |12|. RG 
7| and Monte Carlo simulation IT], [l8|, etc. Among these methods, the RG 



theory has been proved to be very powerful and it has been widely used to investigate the 



2 



critical behavior of different spin systems |2|,l7|,ll9|,|20l. 



flBH 



In this paper, using the RG method, we investigate the quantum spin-1/2 XXZ ferro- 
magnetic model on two kinds of DH lattices with fractal dimensions df = 2.58 and df = 3, 
respectively. Our results show that the systems exhibit finite temperature phase transitions 
in the isotropic Heisenberg limit (A = 0). Besides, we also investigate the above systems in 
the range of — oo < A < and find that they exhibit XY-like fixed points. 

The outline of the remainder of this paper is as follows. In the next Section, the model 
and the calculation method are presented; Sec. Ill gives the results; Sec. IV is a discussion 
about some interesting quantum effects and Sec. V gives a brief conclusion. Some of the 
more tedious of the formulations in Sec. II are illustrated in the Appendix. 



II. MODEL AND CALCULATION METHOD 

The effective Hamiltonian of the XXZ model can be written as 

H = Kj2[a-^) (^r^j + ^M) + <^'^ , (1) 

where K = J/ksT, in which J is the exchange coupling parameter ( J > and J < 
correspond to the ferromagnetic and the antiferromagnetic model, respectively), ks is the 
Boltzmann constant, and T is the absolute temperature, a" {a = x, y, z) are spin Pauli 
operators on the site i. The sum is over all the nearest-neighbor spin pairs and 
A G (— oo, 1] represents the anisotropic parameter. Note that the Hamiltonian contains, as 
particular cases, the Ising model (for A = 1), the isotropic Heisenberg model (for A = 0) 
and the XY model (for A = — oo). 

The DH lattices are all constructed by an iterative manner. Fig. 1(a) illustrates the 
first three construction stages of DH lattice with fractal dimension df = 2.58 (lattice A, 
for simplicity). As can be seen, the initiator is a two-point lattice joined by a single bond 
(construction stage n = 0). Then the initiator is replaced with the generator ( the cluster of 
n = 1 stage). Replacing every single bond on the generator itself, we get the second stage 
of the lattice. If this procedure is repeated infinite times, we then construct a DH lattice 



with self-similar structure |2 ll-l24| . Using the same procedure as that of lattice A, we can 
construct another kind of DH lattice (lattice B, for simplicity) with fractal dimension df = 3 
(see Fig. 1(6)). 



In this section, using the RG method proposed by Caridejl^. [2^, we give the main 
calculation procedure of the XXZ ferromagnetic model on lattice A. The case of lattice B 
can be solved in the same way. 

In the next RG procedure, the Hamiltonians, except for the Ising limit (A = 1) and the 
high temperature limit {K = 0), among the neighboring generators do not commute with 
each other, which leads to the impossible of decoupling the genorator from the whole lattice. 
In order to achieve this goal, the noncommutati vity amon g th e neighboring generators were 



neglected. This method had been used in Refs. [ij, [l5|, |25|, l26| and the quantum effects and 



26| 



the approximation of this method had been detailedly discussed in Ref. 

Based on the above approximation, we take out the generator from lattice A to perform 
the RG transformation, which is shown in Fig. 2. As can be seen, after summation of the 
internal spins a^, a 4 and as (decimation), the generator (Fig. 2(a)) is transformed into a 
new structure, which contains two spins (i.e., (Xi and (72) joined by a single bond (Fig. 2(6)). 
This procedure can be described as 

Tr3,4,5 exp (i?i3452) = exp {H'-^^) , (2) 

where -^13452 and are, respectively, the Hamiltonians associated with the clusters (a) 
and (6) in Fig. 2, Tr3 4 5 denotes the partial trace over states of the internal spins a^, and 
(J5. According to Eq. ([1]), the Hamiltonians -^13452 and are 

i/i3452 = K{l-/\) [{alal + a\al) + + a\al) + (a?a^ + a\al) 

+ {alal + alal) + {ala^ + alal) + (a^a^ + ^a^)] 
+K {alal + a{al + alal + alal + ala^ + alal) (3) 

and 

= K'[{1 - A'){ala^2 + o^l^l) + ^1^2] + ^0, (4) 

respectively, where Kq is a constant included to make Eq. ([2]) possible. 

By calculating the partial trace in Eq. (j2]), we can obtain the RG recurrent relations 
between the new parameters {K', A') and the original parameters {K, A). Firstly, we 
expand exp {H[2) as 

exp (Hl^) = a' + hl^ialal + a^al) + d^^^la^. (5) 
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Note that, in the expansion of exp{H[2), the anticommutation among the spin Pauh oper- 
ators are considered. Further, in the basis which simultaneously diagonalize and (J2, we 
express both sides of Eq. in the form of matrix. The left-hand side of Eq. can be 
expressed as 

/ e^'i \ 



exp {H[2) 



i(e^2+e^3) i(e4_e^3) 



in which 



and 



are eigenvalues of H'12, where 



V 



X' = -K' + 2W' + Kq 



a; = -K' - 2W' + K. 



W = K'(1- A') . 



(6) 



y 



Following the same steps as above, the right-hand side of Eq. 
the form 

/ a' + c'ls 



(7) 
(8) 

(9) 
(10) 

can also be expressed in 



\ 



V 









a — c 



2h' 
12 ^"12 



2b' 



12 



a' -c' 



12 












(11) 



a' + c'ls / 



By combining Eqs. (!6|)- (fTT]) . we get 

a' + c'la = exp {K' + Kq) 



a' — c 



^2 - -(exp {K' - 2K'A' + Ko) + exp {-3K' + 2K'A' + Kq) 



and 



= -(exp {K' - 2K'A' + Kq) - exp {-3K' + 2K'A' + Ko) 



(12) 
(13) 

(14) 



Obviously, the coefficients a', b'12 and c'12 are all functions of K' and A'. 

So from Eqs. ( IT2|) -( !T^ . we obtain the relations between the new parameters {K', A') and 
the expansion coefficients («', &'i25 c'12) follows 



exp {AK') 



(a'-cy^-46;V 



(15) 



exp (4K'A') = + (16) 

Analogously, exp (-^13452) can also be written as 

exp (-^13452) = a 

+ E ^■.Ma.Vla^af, (17) 
(ij(>*)>^(fc,'(>fc)> 

where a, Qj, (ijj,fcz, Cij^ki and /ij^fci are all functions of K and A, which can be determined 
by diagonahzing both sides of Eq. ( IT71) in the same basis of af, a'2, a\, o\ and a\. Because 
the matrices of the two sides are too lengthy, we just give the diagonal matrix form of -^13452 
in the Appendix. 

From Eqs. ([2]), ([5]) and (fT7|) . we can obtain the relations among the expansion coefficients 

a' = 8a, (18) 
&'i2 = 8^12 (19) 

and 

c'i2 = 8ci2. (20) 
These expressions, together with Eqs. ( |T5l) and ( |T6|) . we finally get 



exp(4irO= , '7^' , (21) 



(a + C12) 
la - C12J" - 45^2' 



exp(4A-'A')^ + (22) 

(a - C12 + 2612) 

Since a, 612 and C12 are all functions of K and A, the RG recurrent relations between the 
new parameters [K' ^ A') and the original parameters [K^ A) are determined by Eqs. f l2T|) 
and ( 122|) . However, the analytical expression about their right-hand sides are difficult to 
obtain. In order to obtain the phase diagram and the fixed points, we will numerically solve 
the above two equations and the results will be presented in the next section. 
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III. RESULTS 



By numerically iteration, the RG recurrent equations (12T!) and fl22|) . the phase diagram 
and the fixed points for the system with fractal dimension dj = 2.58 can be obtained, which 
are shown in Fig. 3. The curve (a) in Fig. 3 gives the critical line of the system. As can be 
seen, in the range of A G [0, 1], the phase space is divided into the paramagnetic phase (P) 
and the ferromagnetic phase (F) by the critical line. In this range, we can also find that the 
system exhibits two unstable fixed points, i.e., the Ising unstable fixed point (IP for short) 



[A,l/K) = (1,2.77) and the isotropic Heisenberg unstable fixed point (¥. 



fact, the value of the IP is in accordance with that of the classical model 
the correlation length critical exponent can be calculated by 



21 



P) (0,1.55). In 



2^. At the IP, 



In 6 , , 

= i^. (23) 



in which, b — 2 is the scaling factor and A — i^-^' /^-^)a=i k=0 36 — 

1.86^ 

1.12. We note that this value is the same as that of the classical Ising model 21 



hus, z/i 



23|. It is 



worth mentioning that, at the HP, the system exhibits a finite temperature phase transition 
with critical temperature Tq = ^/Kq = 1.55. However, this is disagree with the result 
of the Wheatstone-bridge-basis hierarchical lattice jl^, where only zero temperature phase 



transition exists, i.e., Tc = when A = (curve (6) in Fig. 3). In addition, the isotropic 
Heisenberg critical exponent is obtained^as = 2.04 which differs greatly from the results 

and other hierarchical lattices with lower fractal 



[z^H = oo) of the 2D regular 



dimensions {df < 2.58) jl4j. Il5|. 



attices 



We also investigate the above system in the range of — oo < A < 0. Our results show 
that the system exhibits an unstable fixed point, which is not given in Refs. 0,0 and we 
call it XY-like fixed point (XYP). The critical line between the HP and the XYP divides 
the phase space into paramagnetic phase (P) and ferromagnetic phase (F) as well. 

Using the same calculation procedure as the lattice A, we can also investigate the critical 
behavior of the XXZ model on the lattice B, i.e., the DH lattice with df = 3. The phase 
diagram and the fixed points are presented in Fig. 4. We can see that the phase space 
is divided into paramagnetic phase (P) and ferromagnetic phase (F) by the critical line. 
In the range of A G [0,1], the system also exhibits two unstable fixed points, i.e., the 
Ising fixed point (IP) and the isotropic Heisenberg fixed point (HP). At the IP (A, 1/K) = 
(1,3.79) we obtain the critical temperature Tc = 3.97 and the critical exponent ui = 0.95. 



Compared with the results of three dimentional Heisenberg model calculated by other RG 
methods^, I27I, o ur results are more consistent with that from series expansion (Tc = 
4.54, z/i = 0.63)|l2j. At the HP (A, l/iT) = (0,2.36), the critical temperature Tc = 2.36 



is very close to that of three-dimensional Heisenberg system {Tc = 2.41) ll|. In addition 



the critical exponent is z/h = 1.69 which agrees well with that obtained by Sousa et al 



on the simple cubic lattice where i^h = 1-64 28|. In summary, we can see that the lattice 



B can be regarded as an approximation for the simple cubic lattice. Besides, just as the 
case of the lattice A, the system on lattice B also exhibits a XY-like fixed point (XYP) 
(A,l//s:) = (-1.99,8.44). 

IV. DISCUSSION 

In this section, we discuss the possible effects of quantum fluctuation. We assume the 
generator of lattice A as a whole. In this case, there will be no noncommutativity and the 
result should be exact in our calculation. With this method, the exact calculation of the 
DH lattice will be noted [K\K , A) , /S! {K , A)]. On the other hand, for the RG transfor- 
mation in Fig. 2, we apply a modifled Migdal-Kadanoff method, which is an approximate 
one. In this method, the original cell (Fig. 2(a)) can be considered as a combination of 3 
arrays in parallel, each of which is made up of two interactions in series. The renormalized 
interaction K and anisotropy A can be flrstly calculated for each combination in series and 



then combined in parallel. The approximate ca^ 
We use the convenient ratios introduced in Ref. 



cu. 



25| 



ation provides [3if (iT, A), A"(ir, A)]. 



and 



^ - K\K,A) ^^^^ 

_ A"iK,A;K,A) 
^ ~ A'(ir,A) ^^^^ 

The T-dependences of and R^ for typical values of A are indicated in Fig. 5. In the high 

temperature limit, both R^ and R^ tend to unity for all values of A ( < A < 1 and A < 

). In the range of low temperature, where the quantum effect tends to drive the system 

to a disordered phase, both R^ and R^ show oscillational behaviours. The fact is usually 

due to quantum fluctuation which, at low temperatures (T) and smaller anistropy (A) are 

important {3] . As we see in Fig. 5, the quantum fluctuation of XY-like model (A < 0) 
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is stronger than the anistropic Heisenberg model (0 < A < 1) at the low-temperature 
region. So, there will be considerable error when we calculate the XY-like model at lower 
temperatures. However, in the range of < A < 1, it is a very good approximation and 
when the anistropy A is negative, as we can see in Fig. 5, it follows the same tendency 
as the positive one. So, we can obtain approximate result at bigger temperatures and give 
qualitatively correct picture of the phase diagram. When it comes to the lattice B, both 
and have the same type of behaviour and we can obtain the same conclusion. 



V. CONCLUSION 



In conclusion, using the RG method, we have studied the quantum spin-1/2 anisotropic 
XXZ model on two kinds of DH lattices with fractal dimensions dj = 2.58 and 3, respectively. 
Phase diagrams, fixed points and critical exponents of the systems are obtained. The results 
show that, in the Ising limit (A = 1), the values of the fixed points and the critical exponents 



agree well with those in Refs. 12|, |22|, |23| . In the isotropic Heisenberg limit (A = 0), there are 



finite temperature phase transitions on the above two lattices. Furthermore, the systems 
exhibit XY-like fixed points in the range of — oo < A < 0. The quantum effects of this 
system show that, at low temperatures, the XY-like model has stronger fluctuation than the 
anistropic Heisenberg model. So, we conclude that there will be a considerable error when 
we calculate the XY-like model at low temperature. Besides, our results also indicate that 
the DH lattice with df = 3 can be regarded as an approximation for the simple cubic lattice. 
As a comment, this method can be extended to investigate spin systems with S* > 1/2 on 
other lattices and we are presently working along these lines. 
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In the basis of al,a2,cr^, cr% and erf, i?i3452 can be expressed as a 32 x 32 matrix in 
form of 



in which A — GX, 



H 



13452 



/aooooo\ 

5 
C 
C 
S 
\OOOOOAy 



c 
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where = X (1 - A) . 
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Figure Captions 



Fig. 1 The first two construction stages of the DH lattices, (a) the DH lattice with 
fractal dimension df — 2.58, (6) the DH lattice with dj — 3. 

Fig. 2 The procedure of the RG transformation. After a step of the RG transformation 
the generator (a) is transformed into a bond (6). 

Fig. 3 Phase diagram of the DH lattice with d^ = 2.58. The critical line (a) separates 
the phase space into paramagnetic phase (P) and ferromagnetic phase (F). The IP, HP and 
XYP, respectively, denote the Ising, isotropic Heisenberg and XY fixed points, {b) is the 
critical fine of the Wheatstone-bridge-basis hierarchical lattice for comparison. 

Fig. 4 Phase diagram of the DH lattice with df = 3. The IP, HP and XYP, respectively, 
denote the Ising, isotropic Heisenberg and XY fixed points. P and F correspond to the 
paramagnetic phase and ferromagnetic phase, respectively. 

Fig. 5 Thermal dependence of the ratios and respectively defined by Eqs. (24) 
and (25), for typical values of A. 
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